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Within the Pauli-Villars regularization technique the fermion propagator is studied in the frame- 
work of Schwinger-Dyson equations in the Euclidean space. Making the generalization of Fukuda and 
Kugo proposals, the analytical continuation is performed into the timelike region of fourmomenta. 
The massive gauge boson is considered and the fermion propagator pole structure is discussed in 
detail. 
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I. INTRODUCTION 
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The natural nonperturbative framework for the study of infrared properties of QCD Green's functions in the whole 
Minkowski space is the formalism of Schwinger-Dyson equations (SDEs). The complex of SDEs is an infinite tower 
of coupled integral equations which must be truncated to be tractable in practice. Most of the studies concerns 
the solution of SDEs in Quantum Chromodynamics pj where it is believed they can give the correct description 
of chiral symmetry breaking and confinement. It is well known fact the the masses of the particles can be largely 
enhanced in strong coupling theory like QCD. The property and difficulty to extract a correct picture of dynamical 
mass generation attracts much attention in recent modeling of QCD Technicolors [H and strong coupling models 

I ' with spinless scalar bosons [5| . The models in two later cases represent the dynamical alternatives to Standard model 
\Q . Higgs mechanism. These models can break symmetry but avoid the introduction of negative scalar boson mass term 
7-H ' for this purpose. 

In QCD, this is certainly the asymptotic freedom and its associated regulating character of the ultraviolet modes 
which makes chiral symmetry breaking physically meaningful. The same is assumed to be valid for Technicolor, 
however where the typical scale A is readjusted about a few orders of magnitudes. In possible strong interacting Higgs 
extensions of Standard model the mechanism can be different since it could follow the ultraviolet completions of such 
' models. In general, the mechanism is not completely understood but it is expected it can be quite exotic as mother 
^H, strong interacting theory can be Q . Having no direct experimental hints which mechanism of mass generations the 
O ' Nature has chosen, we can construct models based on effective interactions. The should describe correctly physics 
O . until some energy scales, above this, the underlying but unknown theory is assumed to be valid. 

■ In this paper we make a look on a simply approximated gauge theory of unspecified symmetry and regularize the 
high energy mode by Pauli-Villars regulators. The Pauli-Villars regulators then can be interpreted as the effect of 
physically hidden sector which effects is observed rather indirectly in the physical spectrum. Very recently, the old 
fashionable idea of Wick and Lee Q has been reopened in Q . To this point we consider strongly interacting theory 
and extend the technique of Fukuda-Kugo solution 0] of gap equation to the case of Pauli-Villars regularization. 
We consider sufficiently heavy gauge boson with mass \x satisfying n < A, where A is Pauli- Villars regulator. To 
see, what is the physical pole mass of fermion which received its mass dynamically, we make a continuation of the 
Schwinger-Dyson equation to the timelike axis of momenta p 2 > 0. We found rather nontrivial results depending on 
the mass ratio /u/A. 

The layout of the paper is as follows. Section II is devoted to the overview of Fukuda-Kugo (timelike) ladder fermion 
SDE with hard cutoff regularization. The results with massless and massive gauge boson propagators inside the loop 
of SDE are obtained by the direct solution of the original integral equation. The absence of fermion poles has been 
actually confirmed in the region predicted by the authors of 0, however, in the light of our further discussion we differ 
in the interpretation. In Section III we do the same for the equation with Pauli-Villars regulator. The numerical 
solution is shown in Section IV. Contrary to the naively regularized ladder QED, the real fermion pole has been 
always identified, at least for nonvanishing (i. Note, the case of vanishing boson mass is more complicated issue and 
to this point the reader can see rather recent paper [Io| where the study of propagator of confined quarks has been 
attempted. Note for completeness and bearing in mind the fact that even gluons can reach their masses dynamically 
[l2| , the technical method proposed here can be straightforwardly used for " QCD with massive gluons" . In Section 
V we discuss the possible form of singularities in the context of Lorentz invariancc constraint. In Conclusion (Section 
VI) the basic results are summarized and further directions of research within this approach are outlined. 

In this paper we will use the following conventions: the positive variables x, y will represent the square of momenta 
in the whole Minkowski space, i.e. for instance x — p 2 for timelike momenta when p 2 > while x — —p 2 for 
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p 2 in spacelike region. Note, our metric is g^ v = diag(l,—l,—l,—l). For purpose of clarity we label the mass 
function B as B s in the spacelike region of fourmomenta and as B t when evaluated for timelikc fourmomentum (i.e. 
B(p 2 ) = B s (x)e(-p 2 ) + B t (x)6{p 2 )). 



II. M ASKAWA-NAKA JIMA /FUKUDA-KUGO EQUATION 

In this Section we review some basic facts on the Maskawa-Nakajima/Fukuda-Kugo equation. This equation rep- 
resents ladder approximation of the fermion SDE with the momentum integration is regulated by an upper boundary 
cutoff A = p% tmax . _ 

In parity conserving theory the fermion propagator S can be characterized by two independent scalar function A, B 
such that S(p) = [A fp — B\~ x (bare fermion propagator is Sq = [/p — m ] _1 ). In the ladder approximation, the 
equation for inverse of S can be written 

£>A{p 2 ) B(p 2 ) =j> m - ig 2 J -0^ la G^(p - q)S(q) 7 p (2.1) 

where mo is a bare mass and G a/3 is boson propagator. The equation Eq. (|2.1j) and the classification of the solutions 
has been discussed in [ll[ especially in Landau-like gauge for which the massive propagator reads 

G af3 (q)= 2 2 J7 . (2-2) 

Few years later, Fukuda and Kugo have found the solution for the timelike momenta for /x = 0. They observe no real 
pole for resulting fermion propagator. Hence the absence of free propagating mode has been interpreted as a sign 
for confinement at that time. First, disregarding the known deficiency of integral cutoff regularization scheme, we 
overview the method of solution and resolve (|2.ip also for nonzero fi. Up to our knowledge such solution has never 
been explicitly published in the literature. We leave the question of the reliability of 'would be' confining solution 
into the discussing section IV.. 

In the next we will follow the paper Q, using A = 1 approximation and making Wick rotation and angular 
integration (see Appendix A) we get 

B s (x) = mo + —7- \ dy — K(x, y, \x ) 



4 Jo y + B 2 (vY 

K(x, y, = *\ , = (2.3) 

x + y + fi z + y/(x + y + [i z Y — Axy 

where A is the cutoff and where C = 3g 2 C2(R) /47r 2 , C2(R) denotes the Casimir invariant of the quark representation 
(C*2(i?) = 4/3 for QCD). The equation (|2.3[) can be quite easily solved by the method of iteration avoiding thus 
inconvenient conversion to the nonlinear differential equation. When mo = but B ^ we talk about dynamical 
chiral symmetry breaking. 

Having kept the solution B s (y) for spacelike y we can define the 'synthetic fermion mass' B{x) at timelike axis of 
fourmomenta such that B(x) = B s (—x) is a solution of integral (|2.3p for timelike x: 



3C f A B s (y) 
4 J y y + B 2 (yY 



B(x)=m + —I dy - | - ^ K(-x,y,fj.) 



K{-x, y, M ) = . 2V 2 ^ 2 — . (2.4) 

~x + y + fi A + y —x + y + + Axy 

which represents correct continuation until the first singularity is crossed on the timelike axis of p 2 . The continuation 
on the usual physical cut of timclike momenta x > (m + fi) 2 reads 

3C f(^~^ 2 Bt(v) 

B t (x) = m Q + B(x) - — dy Jf{, ■ X(x;y,fi 2 ) (2.5) 

4 jo y~ r5 t \y) + l£ 

where X(x; y, /z 2 ) is the discontinuity of K(—x, y, fi 2 ) over the physical cut: 

vi 2n ^(-x-y + fi 2 ) 2 ~4:xy 

x{x,y,v) = - , (2-6) 



FIG. 1: Dynamical mass function for spacelike (left) and timelike (right) region of fourmomenta for massive boson case and 
hard UV cutoff. The solutions are labeled by a = g 2 /4ir. On the right panel the lines representing the imaginary part of B 
becomes nonzero above the threshold. Notice, there are two solutions for the coupling a = 1.4. The pole mass are identified 
by the crossing of thin dotted line (the function y/{p 2 ))- 



and note the presence of Feynman infinitesimal ie, which has been omitted in the original paper 0]. 

To see this continuation was correct let us assume the smoothness of B at the vicinity of the pole y = to 2 
B 2 (m 2 ) — to 2 . Then one can use the functional relation 
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and obtain the dominant contribution to the absorptive part of mass function B 

3C 

ImB(i) = — mirX(x; m 2 ,^ 2 )Q(x — (fi + to) 2 ) . 



(2.7) 



(2.8) 



where stands for the standard Heaviside step function. Clearly, 
when C is a small parameter, 3C/4 << 1. 

In general one have to solve the Eq. (|2.5[) . Assuming smoothness of the function B the equation 
rewritten as follows: 



represents one loop perturbation theory result 

3 can be 
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B t (x) = m + B{x) ~i— ^HjY 
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- 2B t B' t 



X(x; m 2 ,Li 2 )Q(x - (fj, + m 3 ) 2 ) - I(x) 



I(x) 
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dy- 



Bt(y) 



-X{x-y^ 2 )Q(x-^ 2 ) 



4 "Jo " a y-B 2 t { y y 
where j runs over the roots of Eq. y — B 2 {y) — and where we have used the following abbreviation: 



\l-2B t B' t \ 3 = 



l-2B t (y) 



dB t (y) 



dy 
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dB t (y) 



d(yV 2 ) 



(2.9) 
(2.10) 

(2.11) 
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The Eq. (|2.9p represents integral inhomogeneous equations (even for zero Too) with the one singular kernel in the 
integrand of /. It can be solved by a standard numerical method and we provide some details in the Appendix B. The 
dynamical chiral symmetry breaking is of great interest for us and we will concern on these solutions in this paper. 
Inclusion of small explicit chiral breaking term is straightforward and we leave this case aside of our interest. 

In order to scale the dimensionfull quantities we take fi — 1 in arbitrary units. Further, the numerical hard cutoff 
is adjusted to be A = e 5 /i ~ 148/1 . Its value constraints the meaningful region of the square of fourmomenta as 
p 2 G (—A 2 , A 2 ). The numerical solutions for various coupling constant are presented in Fig. 1. 

As expected already from the discussion in Q there are several phases of the theory. They have been confirmed 
and actually found also here. They are as follows: 
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I. Chiral symmetric phase- for the coupling bellow the critical value a c , there is only trivial solution B = 0. 

II. Chiral symmetry breaking phase, where the propagator has a real pole. In this case we get also nontrivial 
solution for B. For our parameters fx, A this phase is characterized by the coupling strength a £ (a c , 1.4) and the pole 
mass m is at most of order fj,. In fact, we can observe and have found two solution for a = 1.4. Although numerically 
they both cut the axis y/p 2 , the one of them only sweep the region under the diagonal and their timelike high energy 
asymptotic crucially differ. Of course, the question of uniqueness of such continuation naturally arises. 

III. Chiral symmetry breaking phase above a — 1.4, where the mass ratio m/fJ, ~ 4.2 (and B(0) ~ 3/i) is achieved, 
the pole and associate analytical cut vanish and the mass function blows to 'infinity' like \f[p 2 ) for large p 2 . Note, 
no sign of such transition is observed in the spacelike regime. Clearly the spacelike solution does not reflect dramatic 
changes that happen in timelike part of the Minkowski space. The solutions in the regime III. have been originally 
called 'confining' because of the absence of real pole and associate particle production threshold. In fact we have to 
be more careful with the interpretation, particularly because the Wick rotation is not allowed, the rotation of the 
contour cross the singularity at timelike momentum infinity. Clearly, in this case the hard cutoff regularization of 
the momentum integral does not commute with the Wick rotation and there is neither confidence that we obtain a 
true solution in the whole Minkowski space. This is because we are dealing with non-asymptotic free theory where 
the result are not independent on the UV integral cutoff. We argue here, the "would be" confining solution would be 
meaningful only when the theory is defined in the Euclidean space from the beginning. Clearly, this is not the case 
of ladder QED and the observed nontrivial solution can be an artifact of given scheme. Although we suppose that 
these things are rather well known, we discussed them explicitly here for the purpose of clarity. 



A. Massless vector boson 



The solution for massless vector boson has been firstly obtain in Q. For this purpose the authors of the paper 

Q transform the integral equation into the differential one which has been solved numerically. Here we confirm 

their solution by direct solution of the original integral equation. Taking a limit fi — ► in Eq's (|2.3[) and (12. 9|) is 
straightforward. The SDE for selfenergy in the massless case (mo = fi = 0) then reads 

n ( \ a (V A Bs ^ , < >A , 010 s 

B s (x)= dy[--l) + (2.12) 

Jo v ^ / y + Bj(y) N c 



B t {x) = —I dy[l--\ 2.13) 

N c 4 J V xJ y-Bf(y) 

where N c is the number of colors and where we have already omitted intensifimal imaginary prefactor ie since as it 
is known that the equation y — B 2 (y) = has no real roots in hard cutoff theory. The inhomogeneous term in (|2.13p 
is simple constant-the usual fermion condensate: 

b = <^ >a = ^ r 2 (2 . 14) 

N c 4 J y + B^y) V ' 

The numerical results are shown in Fig. [5] for three distinct values of the coupling constant. Similarly to previously 
discussed region III of massive boson case the results could be interpreted with a great care. Again the Wick rotation 
invalidates because of naive regularization scheme and the resulting timelike high momentum behaviour can be an 
artifact of inappropriate calculation scheme. 



III. DYNAMICAL MASS GENERATION WITH PAULI-VILLARS REGULARIZATION 



In this Section we consider the fermion gap equation regularized by Pauli-Villars regulator function. Such construc- 
tion is very straightforward in our approximation since it achieved by the replacement of the gauge boson propagator 
by the sum of ordinary , i.e. original one propagator, and the "wrong sign" boson propagator but with Pauli-Villars 
mass A instead of /i. Explicitly, we make the following in the kernel of the gap equation 
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(3.1) 
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FIG. 2: Fermion dynamical mass function in ladder ap- x 
proximation with hard cutoff A 2 = e 10 (in arbitrary 
units) for the three value of the coupling constant a = 

1.4, 1.5, 10. The thin dotted line represents \ft?, the so- FIG. 3: Dynamical quark mass for a = 4 as described in 
lution for timelike momenta are labeled by thick 'upper' the text, 
lines, the spacelike solutions are added for the compari- 
son (decreasing lines of the same type belongs to the same 
coupling, clearly B s and Bt are identical for x = p 2 — 0) 



Within a finite A we do not consider the hard cutoff, which is formally removed from now. 
and where we have also introduced effective gluon mass fj,. 

In order to safely define the pole fermion mass we consider massive boson case in presented study. We suppose 
that the timelike singularity can be drastically changed due to the presence of strong interaction and we leave this 
complicated problem for separate numerical search and discussion elsewhere. However, as usually we assume no 
dynamical singularities in the fermion propagator in the first and the third quadrant of complex p 2 plane. It not even 
allows to perform standard Wick rotation, but the main advantage of presented method is that we can readily follow 
the Fukuda-Kugo trick and make the continuation to the timelike axis. 

Doing this explicitly we get for spacelike solution: 

B s (x)=m + - I™ dy fffi K(x,y) 
K{x,y)=[K{x,y^ 2 )-K{x,y,A 2 )\ 

= ^-{n 2 -A 2 - y/(x + y + ix 2 ) 2 - 4xy + y/(x + y + A 2 ) 2 - 4xy} (3.2) 
Ix 



and the continuation on the analyticity cut reads 

t \ t>, s Tt \ , ■ ^Xix-^m^^eix-ip + m^-Xix^lA^Qix-iA + mj) 2 ) 
B t (x) — B(x) - I(x) + ta x - — : 
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dB t (y) 



(3.3) 



vr 



B(x) =m + - f°° dy f'ff K(-x, y) , 



I{x) = -P. / dy gf^(*; V, » 2 )S(x - fi 2 ) - —P. / dy ^^(x; y, A 2 )Q(x - A 2 ) 

it Jo y-Bfiy) tt J y-Bf(y) 



K(-x, y) = K(-x, y, fj 2 ) - K(-x, y, A 2 ) 
1 



, {[i 2 - A 2 - y/(-x + y + n 2 ) 2 +4:xy + y/(-x + y + A 2 ) 2 + 4xy} . 
Ix 
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FIG. 4: Dynamical fermion mass for a — 5 
and fi/A = v / 5X 



FIG. 5: Dynamical masses for 'super-strong' coupling 
as described in the text. Solid lines represent result for 
a — 10 and fi/A = vOTT, dashed line stands for a = 5 
and fi/A = VttooT 



(3.4) 



The equations (|3.3|) have been solved numerically by the method of iterations, some useful details concerning the 
numerical treatment can be found in the Appendix. 

Remind here known feature of QCD scaling: the constituent quark masses are of the same size as the QCD scale, 
i.e. -B(O) ~ A. In this paper, for any reason, we freely focus on a larger region of parameter space providing us 
the results for softer coupling B(0) < A and stronger B(0) > A couplings as well as. We plot sample of numerical 
solutions in Fig. 3-6. For most of the solutions we fix the mass ratio to be A 2 //! 2 = 10, the exception is explicitly 
mentioned. The critical value of the coupling has been identified a c ~ 3.8, bellow that we do observe the trivial 
solution only. In Fig. 3 the numerical solution is plotted for the coupling strength a = 4.0. Being rather close to the 
critical value a c , the quark propagator has a one real pole at some point m which value is very closed to the infrared 
mass B(0). Without any doubt, the constituent mass can be freely identified with the infrared mass or with the pole 
one. In Fig 4. the solution for a = 5 is shown and we observe that propagator develops two real poles under the first 
branch point. The absorptive part of B is largely enhanced because the both differentiations dB(m)/dm are not so 
far from 1 (see Eq. (|3.3[) ). Such solution is possible since the boson is massive, however it may be an artifact of our 
approximation. The observation of two pole solutions is in agreement with similar observation mad in [l3|, however 
we regard this as a curiosity rather then possible physical scenario. Increasing the coupling further, then the second 
pole vanishes at some point and we retain with the one real pole solution again. This situation is exhibited in Fig 5. 
. Note, the propagator becomes largely enhance bellow m, since the function B t lies very closed to the diagonal of 



■s/ffi, Bt graph. Contrary to the case of solutions with two poles, the imaginary part of B appears to be suppressed 

now, because of differentiation B'(m) < 0. In that case the pole mass largely differ from the infrared value B(0), it 
is typically few times higher. 




IV. DISCUSSION OF THE RESULTS 



In this section we continue the discussion of the results and qualitatively compare with the already known results 
presented in literature. We also discussed the location of the poles in the complex plane and its relations with the 
violation of Lorentz invariance. 
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FIG. 6: Real part of the dynamical mass function 
for x = p 2 . The lines represent the models 

characterized by the couple of parameters a and the ratio ^i 2 /A 2 as shown in the box. 



In QCD with chiral quarks the asymptotic freedom ensures the finiteness of dynamically generated masses without 
any use of regularization scheme. In our case the infinity is regulated by the Pauli-Villars term induced by auxiliary 
unphysical field added to the Lagrangian, hence the mass function would be finite even when the bare mass is added. 
As we have checked explicitly for small bare mass, it affect the ultraviolet behaviour of the mass function but has 
a little effect in the infrared and the singularity structure remains the same. The main result of our paper is the 
numerical observation of the fact that the real singularity is never absent in our translation invariance preserving 
calculation scheme. The singularity is a real pole for large parameter space of a,fi and A However approaching 
the limit of vanishing boson mass fi, the pole is going to meet the branch point and the situation is getting more 
complicated. In Fig. 5 the result is plotted for a small ratio /z 2 /A 2 = 0.001. We should note that the exact massless 
limit may not be well defined in approach presented here so that Fukuda-Kugo continuation method is not fully 
justified in this case. Recall, the functional used in our equations is defined on the set of differentiable functions, while 
in contradiction, the differentiation B'(m) may not exist in the exact massless case [i = 0. The zero boson mass case 
will be studied elsewhere due to its own serious physical consequences. 

Stating without proof, we suppose that true mass generation and the evidence of nontrivial pole of the fermion 
propagator go hand by hand with the property of asymptotic freedom (decreasing of effective coupling at ultraviolet), 
whilst the absence of real pole as happened to the propagator function in Fukuda-Kugo equation is artifact of im- 
proper calculation scheme (hard cutoff regularization). Further note, the numerical results obtained here qualitatively 
correspond with the results of the papers [l3l |. where the same model was studied in the Minkowski space for the 
first time. At this place we discuss some common and some distinct features of solutions presented here and in the 
paper 13]. In both approaches the validity of Wick rotation is assumed. In other words, no (dynamically generated) 
complex singularities are assumed in the first and the third quadrant of p 2 complex plane. The main distinction is the 
analyticity assumption explicitly used in the paper [l3|, where it was assumed that the propagator is holomorphic in 
the whole p 2 complex plane up to a real positive semi-axis and that the propagator satisfy certain integral (generalized 
spectral) representation. To this point, there is no similar assumption in the approach developed by Fukuda-Kugo 
and followed in presented paper. In the method presented here the correlator at timelike regime of fourmomenta 
is built on a base of the knowledge of its spacelike counter-partner. However, being possibly out of the domain of 
analyticity we should be aware that such continuation is by no mean unique. On the other hand , there is also 
reasonable quantitative agreement with [l3|, especially when we are not far form the critical coupling characterizing 
dynamical chiral symmetry breaking. Actually, depending on the coupling strength a the propagator develops one or 
two real poles with corresponding branch points providing their values are in an approximate agreement with [131 ] . 

In the literature there is a certain effort to guess the complex structure of QCD Green's function by the method of 
continuation of Euclidean result or by indirect reading from the phenomenological consequences which would follow 
form particularly assumed singularity structure of Greens functions. Complex conjugated singularities of quark and 
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fermion propagator in planar strong coupling QED have been considered in the context of confinement 14J and in PT 
symmetric Quantum Field theory [15|]. Additionally recent studies have modeled Euclidean space lattice data with 
propagators that have complex conjugate singularities [l(| H3] • More phenomenologically, the meson bound states [H[ 
and the parton distributions (l9j has been calculated with a quark consisting of pairs of complex four momenta. In 
many, if not at all of these recent studies, the fact that complex singularity structure can affect the Lorentz invariance 
of the theory has been overlooked. Therefore, in what follow we mention the question of Wick rotation, location of 
complex poles and possible lost of Lorentz invariance in order to pay attention for. 

Actually, increasing the coupling strength one can expect that apart of the real pole a new complex singularities 
appear. The area of real part of p 2 where we can expect new complex singularities in the propagator function is 
indicated by the ellipses in Fig. O In recent, we are not able to estimate the characteristic of the behaviour and/or 
the position of these complex singularities with good confidence. Instead of this, we would like to present the argument 
which largely enforces our believe in the validity of Wick rotation in a case of full realistic solution of Schwinger-Dyson 
equations. Recall at this place, there is known relation between location of complex singularities of Greens functions 
and the Lorentz invariance of the theory. Indeed, its more than 30 years known that complex singularities located 
simultaneously on both side of the real axis of p 2 automatically generates Lorentz violating peaces in the propagator 
itself [20j . In this case the naive Wick rotation invalidates and one has to account the complex singularities in a suitable 
manner. Actually, assuming the complex poles (with non vanishing Im part) one can exhibit the existence of Lorentz 
violating pieces originally not expected in. The evaluation of the appropriate Feynman integral has been exhibited 
for the case of complex conjugated poles in 20] (some prefactor has been corrected in the paper [2l[, the actual 
evaluation has been performed for scalar field content only, the extension to the loop integral with internal fermion 
links is rather straightforward) due to the rather different reasons. The generalization to more general locations of 
singularities were already discussed in [20| . From these arguments it follows: if the true exact solution respect the 
Lorentz invariance of the theory then selfconsistence of the Schwinger-Dyson equation solution for quark propagator 
can not involve complex poles on the both sides of the real axis of the square of the fourmomentum. In Lorentz 
invariant theory, the Wick rotation is indirectly justified in this way. 



V. CONCLUSION 



The model we have described , is of course not true QCD and hence wc do not claim that the fermion propagator 
we have studied here is truly representative of the singularity structure of the quark propagator. In a more realistic 
models, which include the dressed vertex unction and/or higher order skeleton graphs, one expects that the details 
of the propagator structure would be different. It should also incorporate the aspect of asymptotic freedom in a 
more proper way. However, our study demonstrates that singularity structure of the quark propagator is very likely 
dominated by the real singularity and we expect a quantitative but not a qualitative changes when the approximations 
improve. We discuss possible singularity structure in the entire complex plane of momenta. We argue, in the Lorentz 
covariant formulation of QCD (and any QFT) the singularities of the propagators (and of integral kernel of gap 
equations at all) must appear only on the one side -upper or down- separately for the left and the right half complex 
plane. Otherwise we necessarily sacrify Lorentz invariance of the theory. This argument strongly support the validity 
of Wick rotations in Lorentz invariant theory. 
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APPENDIX A: ANGULAR INTEGRATION 

The following integral: 

f dz ^ l ~ z2 = tt(-o + \Jo? - 1) , (Al) 
J-x z-a 

is used in order to perform the angular integration in SDEs. 
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APPENDIX B: PRINCIPLE VALUE INTEGRATION 



1. / for Fukuda-Kugo equation 



SDE for timelike momentum (|2.9| represents the complex inhomogeneous integral equation with singular kernel. 
In this appendix we describe some details how to numerically deal with. The integral to be evaluated reads 



d y '^r-Xix-v,^)^-^) 



(Bl) 



noting that the other terms in the SDE represent regular integrals with a smooth kernel. 

Let us assume that we have made a good guess of the value of the pole mass to, then it is convenient to write down 
the expression (|B1[) separately for various regime of momenta term / equivalently as follows 



I(x) = 



x < fl 



(B2) 



The function B(y) is real bellow the thresholds and the kernel is regular for x bellow the threshold, hence there is 
no need to denote P. in front of integral since y < m 2 for \f{x) < to + fi. At the point y/\x) = m + fi the kernel 
singularity is suppressed by vanishing function X, thus we can safely write 



I(x) 



3C 
4 



y-Bf{y) 



jtx < x < (/i + m) z 



(B3) 



For a larger x we necessarily cross the singularity in the kernel, however B(y) remains real to the threshold and we 
consider this regime separately: 



I(x) 



3C 



P. 



Ml-tf. B t (y) 

dy nTr^X(x; y, [i , 



(/i + m) 2 < x < (2/i + to) 2 



(B4) 



Increasing x the kernel become complex and we divide the integral to the P. value integration over the real B and to 
the regular integration over the regular kernel with complex B: 



I(x) 
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T 
T 

for 



P. 



(m+M)2 , B t (y) 

dy- — Twrx x {x;y,fj, ) 



y-B 2 ( y y 

dyX(x;y,fi 2 ) 



y - 



Y 2 - V 2 



(y 



x 2 R 



4S|S 2 



i S; 



(y 



4S|S 2 J 



x > (2/i + m) 2 , 
where we have used following shorthand notation 

E fl = ReB t (y); £/ =ImB t (y). 



(B5) 



(B6) 



To avoid some unwanted numerical fluctuations which usually stem from asymmetric distribution of mesh points 
when P. integration is numerically performed, we use a standard trick. Consider for this purpose the third of considered 
momentum regime where x runs over the interval ((/i + to) 2 , (2/i + m) 2 ) (see (|B4jl ) . In our numerical treatment the 
integral is replaced by a discrete sum with the appropriate (Gaussian) weights, i.e. 



I(x) 



3C 



P. 



dy , X{x-y^ 2 ) 



3C 



BtM 



X{x;y v ^ 2 ) 



(B7) 



The numerical fluctuations are subtracted by the following trick 



3Cv- . , B t ( yj ) 



-T^2 w (yj) 



yj -Bfivj) 



X{x\yj,n 2 



BtiVj 



yj-B 2 (y 



_», , sX{x; y:j ,fi 2 ) 



mX(x; to 2 , /i) 
yj - to 2 



+ X{x-m 2 ^ 2 )\og\ 



{■sfx — /i) 2 — TO 2 



(B8) 
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where the later two terms vanishes in the exact continuum limit. We have found this approach is very stable and 
the pole mass m can be identified after few runs of the iteration program. In order to achieve better stability, the 
actual search of the root of equation B(x) — x = has been performed by hand for each solution. After making a 
few iterations in rrij then the achieved numerical accuracy of our search is estimated by Am ~ step/m where 'step' 
means the difference of two neighboring points at vicinity of m. Likewise in the case of SDE in Euclidean space, the 
integrals were cut by the ultraviolet cutoff. 



2. Kernel with Pauli-Villars propagator added 

Generalizing to the case of Eq. (|3.3[) is rather straightforward. The principal value integration is treated in the 
same fashion as in the case of Fukuda-Kugo equation. Now the relevant term reads 

I{x)=I^{x)-I K {x), (B9) 

where 1^ and I\ are the integrals considered previously, but where \x is replaced by QCD scale Aqcd in the later case. 
Contrary to the previous case the integral is insensitive to the value of upper boundary Ah, when it Ah > Aqcd 
For completeness we list the function Ia bellow 



Ia{x) 
Ia(x) 
Ia(x) 







x < A 2 
3C /-(^- A > 2 B t (y) 

(™+^ 2 B t (y) 



4 



A 2 < x < (A + mf 

(A + mf < x < (n + A + mf 



3C 

T 
T 

for 



P. 



(v^-A) 2 



y-Bf( y y 

dyX(x;y,A 2 



i E; 



x> (/i + A + mf 

The treatment of principal value integrals is the same as in the previous case of Fukuda-Kugo equation. 
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